On the mutual polarization of atoms He 4 and a possible nature of 

the electric activity of He II 



We propose a simple method based on the standard quantum-mechanical perturbation 
theory to calculate the mutual polarization of two atoms He A . A possibility of the explanation 
of the electric signals observed in several experiments with helium-II is discussed. 
PACS: 67.40.Pm, 32.10.Dk 
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1 Introduction 

As usual, the superfluid He A draws attention in view of its superfluid (SF) properties. 
However, the experiments [I] [3] revealed electric properties of He II in the cases where 
external electromagnetic fields are absent. These results were rather unexpected in that the 
free He A atoms possess no charge and dipole or higher multipole moments. 

Let us dwell in more details on experiment |T], where a second-sound standing half- wave 
was generated in He II placed in a metal resonator. The electric signal (potential difference) 
oscillating with the frequency of the second sound and the amplitude aU ~ ksAT '/2e was 
observed between the electrode on the resonator wall and ground. Here, e is the electron 
charge, and aT is the amplitude of temperature oscillations in the second-sound wave. 

Although the attempts to understand a nature of this signal were made in several papers 
[4]-[7], its satisfactory model has not been constructed yet. In particular, it was assumed in 
[4] that the quadrupole moment generates a dipole moment (DM) of atoms He II in some way. 
However, neither a mechanism of the DM generation nor a nature of the quadrupole moment 
were explained. In work [5], an idea was proposed that atoms He II acquire a DM due to 
their acceleration in the second-sound wave. However, the amplitude aU obtained in such an 
approach turned out to be two or three orders of magnitude smaller than the observed value. 
In [6], it was suggested that the results found in [TJ are based on the acceleration-induced 
polarization of atoms [5] and the presence of a large enough entropy s s ~ 2ks (per atom) in 
the SF component of helium. A possible explanation of the origin of this entropy is given in 
[7]. However, the initial equations were not quite rigorously deduced in Ref. (6j, so that the 
relation for the second sound 



for s s ~ 2k B and T = 1 2 K yields the value of u 2 by one or two orders of magnitude larger 
than the observed one (although, for s s = 0, this relation agrees with the experimental value 
very well). 

In those works, the main mechanism of the polarization of atoms was the inertial 
mechanism [5]: a polarization due to the inertial forces at the acceleration. However, the 
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inertial forces (and the gravity force) give, as known, the identical acceleration to bodies with 
different masses. Therefore, in our opinion, the acceleration by itself cannot be a primary 
source of the polarization of atoms He (see details in [§]). Works [9J [10] advance the idea of 
the polarization of atoms due to the interaction with one another ("tidal" mechanism). In 
this case, the acceleration can give a correction to the polarization, but in the frame of the 
tidal mechanism. It was assumed in [10] that the polarization of He-II can be related to the 
ordering in an ensemble of vortex rings. However, the value of the DM of a ring was postulated, 
and, in addition, it was found that the volume polarization of any nature cannot explain the 
properties of a signal [I], since a volumetric signal must strongly depend on the resonator 
sizes, whereas the experimental signal is independent of those sizes. The tidal mechanism 
ensures [9] [II], [8] the polarization of helium in the presence of the gradient of density or 
temperature. Such gradients allow a second-sound wave to induce a volume polarization in 
helium and a weak volumetric signal less intense than the experimental one by 1-2 orders [8]. 
Works [12] and [1.3] presented the idea of the relationship of the experimental electric signal 
to the unidirectional polarization of atoms He 4 located at the electrode surface. This idea 
yields an approximate microscopic model [14] which explains significantly well the properties 
of the observed signal. 

The application of the tidal mechanism requires to know the mutual polarization of two 
He 4 atoms. It was determined in [151 HE], where it was also shown that each of two interacting 
atoms He 4 induces a dipole moment 



on the neighboring atom. Here, n = R/R is the unit vector toward the second atom, and 
as = h 2 /me 2 = 0.529 A. The calculation in [T5l EE] was executed on the basis of special 
dispersion relations known mainly to experts in atomic physics. Therefore, we propose to 
solve the problem by a simpler method which does not leave the frame of the ordinary course 
of quantum mechanics and, in this case, will verify the result (DQ). Such a calculation is also 
actual because, in order to fit the experiment [I], it is important to know the tidal DM (pQ) 
with maximum accuracy |14| . 

2 Mutual polarization of two atoms He 4 

According to the standard perturbation theory (PT), we will search for the ground-state wave 
function for two immovable interacting atoms He 4 , A and B, as the expansion in all possible 
states of two free atoms A and B: 



where rriA and m# are the complete collections of quantum numbers (n, /, m) for these atoms. 
By = we denote the set of states with a given I and all possible corresponding m and 
n (m — 0,±1,... ± /, n — l,...,oo). In what follows, we will omit sometimes the tilde for 
simplicity: tyf = \l) = \l), C55 = Coo, etc. Since the excited s-states (lsns) give the zero 
contribution to a DM (Tj]), we neglect them, so that ^0 means the ground state lsls of an 
atom He 4 . It is obvious that = Cji. According to PT, we have 
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where the upper index shows the order of PT. The tidal DM of atom A is sought by the 
general formula 

d A = J ^AB^AB dT A dT A dT B dT B } (4) 

where are the coordinates of electrons, and 

d A = e(rf + 4). (5) 

The function describes a parastate. Therefore, the coefficients are nonzero only for 
parastates tyj. In view of this fact, the spin function is the same for all SSfj in ((2j). Therefore, 
the spins can be took out in (j4j). Since the sum over spins gives 1, we omit the spin coordinates. 

We note that the quantity d A ~ cio(l) ~ (as/R) 4 was determined in [TpJ. The coefficient 
Cio was found in the first order of PT, but its calculation includes an inaccuracy: in the 
charge operator Q (see formulas (11) (13) from [10J), the Z-axis was chosen prior to the 
differentiation. But this axis should be fixed after all differentiations. In this last case, all the 
coefficients c[J become zero, because ~ (0 jB |<5,b|0' b )-R~ 1 and 

(0 B |Q B |0 B )i?- 1 = 0. (6) 
In particular, (0 B \QB, q \0 B )R~ 1 ~ Yl ^ af~ dl~ 1^ = = 0? an d we S e t zero only for the 

a/3 

Coulomb potential. Relation ((6)) is valid, because an helium atom has no moment in the 
ground state. 

In this connection, we will determine C\o in the second order of PT: 
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where m = (ra^^B), and 

(minima) = J ^^U^gndrt dv^dvf drf . (8) 

It follows from U = Q\Q B {1/R) \T7\ and © that (00|f>|00) = 0. Therefore, only the second 
term on the right-hand side in ((7j) is of importance. The main contribution to this term is 
given by the states fh = (1, 2) and m = (2, 1). We have 

(2) ^ (m|fr|00)(10|fr|m) 

m^(0,0) 

As ni and n 2 increase, the terms in j9j) decrease. The calculation indicates that d A receives 
a significant contribution from a lot of states at Cio with n > 2 (in particular, n = 5 and 
6). Since the number of states is huge, it seems impossible to find and to sum over all the 
required states in ((31), as well as all cio with different n and m. 

Nevertheless, this difficulty can be avoided if the series is summed approximately. We note 
that the difference of energies for an atom He 4 is minimum for those of the first excited 
level ls2p and the ground state, aEi{ji = 2) = 21.236 eV, and the maximum difference is 
represented by the ionization energy AE ion = 24.58 eV [18]. In ([9]), a significant contribution is 
given, for example, by the state ls5p c aEi(u = 5) = 24.066 eV [T8] . Since aE are sufficiently 
close for different n, we set them to be equal to a constant aE. 



Let us supplement © by the term m = (0,0) <00|f7|00> (10|£>|00)/2a£ 2 (it is equal to zero 

,(o 
o 



due to (00 \U 1 00) = 0). With regard for the fact that AE lfi = E^ 0) -E[ o) o = 2E^ 0) - E[ 0) - E^ 



E^ - e{ 0) « -aE and aE^ = E$ - E$ = 2Ef ] - E*^ - E^ « -2aE in Q, we obtain 
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where q depends on n and m in the state |1). By writing (fTOl) , we used the well-known 
property of completeness 

(frA = 1, (11) 

m 

where the sum is taken over all states of both atoms. By considering only the terms with coo, 
Cio, and c i in ([2]), relation (J4]) yields the first part of the DM of atom A: 

d A = J drfd^drfdr^dXo^ + + c o?*o*f]| 2 (12) 
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« <So£c$<00|d l 10 ) + c - c - = ^E( °l a |10)(10|t/ 2 |00) + c.c. 
1 1 

Here, ^2 is the sum over all states |1), i.e., over m = 0; ±1 and n = 2, 3, . . . , oo. Besides |10), 

i 

we include all the remaining states (m^me) in (fT2l ): = |lsnp), and also |lsns), \lsnd), 
\lsnf), etc.; and, in addition to |0 B ), all higher states |m). This is possible to make because 

(00 1 d \ftiAfriB) = for them. Then, with regard for (CCD), we obtain 

d? « ^2 E (00|d A |m A m B )(m A m B |f/2|oo) + c.c. = £|L( 0|d A f/ 2 |00>. (13) 

" m A m B 

Relation (fill and the normalization conditions for ^ AB yield 

c 2 -l- V c 2 ~, ! (001^100) _ ! , E Vd y 

Ziji^O.O) 

where E VdV ~ — ^f^^ ~ — 2.15e 2 (a^/i? 6 ) is the van der Waals interaction energy for two 

atoms He 4 . We have |£V<ftH <S therefore, coo ~ 1. 

With regard for (j2J) and (SJ), the total DM of atom A is given by the formula 

d A = £ 4^mf c m^mf (m^mf |d A |m^mf). (15) 



The above-derived expression for di is a part of the last relation from the terms CqqC\q in 
the second order of PT. The formula for d A includes also the nonzero contribution in the 
first order of PT from terms of the form c^c^ , c^j , and higher-order ones which behave 
themselves, respectively, as ~ (ae/i?) 7 , ~ (a^/i?) 9 , etc. (see below). At I and j which are not 
equal to zero (0) simultaneously, the relation 

(1) _ (l A j B \U\0 A B ) 

j E^ ] - E {0) + Ej 0) - e\ 0) 



is valid. Since (lj\U\00) ^ only at I, j > 1, we set 2Ef ] - E)"> - E- t 
contribution in the DM from c 0Q is taken into account in di, we redefine c o in (fT5l) according 
to f fT6l) with —2aE in the denominator (such c 00 is equal to zero). Since c^q = 0, we will write 
also —2aE in the denominator for it instead of —aE. We can now use the relation (fTTT) . and 
relation ffT5l) minus di is reduced to 
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Whence we obtain the total DM of atom A in the principal (~ (a^/i?) 7 ) approximation: 



d A = d? + d£ « 



(0 A K |d'[/ 2 |0 A s ). 
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Using the formulas for the operator of charge Q [10], we represent the interatomic potential 
U as the sum 

U = QaQb^ = Udd + u qd + u qq + ... 

where 



U. 



del 
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(19) 
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is the dipole-dipole part of the potential, and 
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+ 4d B nX:Qi-(A^fi) 



(21) 



is the dipole-quadrupole part. The quadrupole-quadrupole part, U qq , includes a negligible 
correction to the DM (~ (clb/R) 9 from U qd U qq ). Therefore, 



Uj d + 2U dd U ad + Ul, 



dd^gd 



q d' 



(22) 



In view of the parity relative to the inversion, we preserve only the term 2U dd U qd . As a result, 
we obtain 

K A 

(23) 



(o A o B \d-u dd u qd \o A o«). 
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Let us choose the Z-axis along R: R = Ri z . Then a nonzero contribution to d is given by 
the terms 

d% d U qd = ^ {{d B x fd A d A Qt + (x~y)- {d B )\d A f{Q A x + Q A 

With regard for the isotropy \I/o = |0), we eventually arrive at the general formula 
30L 



(24) 



(o B \(d?f\o B ) ( (o A \d A d A Q A z \o A ) + (o A \(dmQt - QtW) 
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which can be rewritten in the following simple form: 



-D 7 e 



(25) 



(26) 



Dr = Z A Z B (^j (ri) B (rt)A. (27) 

Here, (ff) B = (0 B |(rf /a B ) 2 |0 B >, (rf) A = (0 A |(rf /a B ) 4 |0 A >, Ry = e 2 /2a B = 13.6 eV. 

The next correction to d A is of the order of ~ (a B /R) 9 , and it is ~ (R/a B ) 2 times less 
than (l26l ). For He II, i? ~ 3.6 A, (R/a B ) 2 — 50; therefore, it can be neglected. It is of interest 
that the separate corrections to D 7 for various n and m in the (11) and (12) states from (j2J) 
are small, about 0.1 and less. Moreover, we have an alternating series for cU, so that several 
hundreds of various corrections and about 10 different n significantly contribute to the total 
DM D 7 (27)). 

The positive sign of D 7 means that the electron cloud of each of two atoms He 4 shifts 
to the neighboring atom. That is, the interatomic space acquires an excess of the negative 
charge. 

To calculate the value of D 7l we need to know the wave function |0) of the ground state 
of an atom He 4 . The exact analytic solution is not available for it, but some approximations 
have been proposed. Within the simplest one-parameter one [321 CES] 



2 

xx2\ /~4 



= ^(ri)^i.(r 2 ), <p ls (r) = Z*yna 3 B e' rZ ^ aB (28) 

(Z* = Z - 5/16), we get (f 2 ) = 3/(Z*) 2 « 1.053, (f-f) = 45/2(Z*) 4 w 2.775. In ([271) for aE, 
we take the least (aEi(ti = 2)) and largest (AE ion ) values, so we obtain D 7 « 14.32 4- 19.18 
(here we take into account also that Z A = Z B = 2). 

The more exact multiparameter \&o were calculated within various methods. In [20], it was 
proposed to seek in the form 

oo 

V (s,t,u) = (p(ks,kt,ku), <p(s,t,u) = Ae' s/2 ^ c n , 2 i, m s n t 2l u m , (29) 

l,n,m=0 

where 

s = T\ + r 2 , t = n — r 2 , it = r 12 . (30) 
The formalism in the variables s, u, t is given in [21] in detail. For the 3-parameter \l/ [20l l2T] . 

up = Ae- s/2 (l + cm + c 2 t 2 ), ci « 0.081, c 2 « 0.01, jfe « 3.63, (31) 

the calculation yields 

A = (A; 3 /4vr)(4 + 35ci + 48c 2 + 96c 2 + 576c 2 + 308cic 2 )" 1/2 « 1.325, (32) 

(f 2 ) = 768 ^ 2 (l + + 24c 2 + 42c 2 + 480c 2 + 195 Cl c 2 ) , (33) 

. 4 , 23040tt 2 A 2 / 979 r „ „ n , 1W „ . 8305 \ , . 

( f l) = 7J0 ( 1 + -^ C l + 52 °2 + 68C ? + 1704C 2 + ^ ClC2 ) • ( 34 ) 

We have also performed the numerical calculations for 6-, 10-, and 20-parameter Hylleraas 
functions [201 [23] and for more general 10-, 39-, and 80-parameter Kinoshita functions 
[22 [25] of the form 

oo 

V (s,t,u) = <p(ks,kt,ku), <p(s,t,u) = Ae' s/2 ^ c n ^is n+l (t/u) 2 \u/ s) m (35) 

l,n,m=0 

which include the negative degrees in the expansion. The results are given in the Table. For 
the 3-parameter \l/o, the result is obtained numerically and analytically. It is seen from the 
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Table that the more the number of parameters, the less the variation of the results for (ff), 
(ff), and D 7 . Moreover, the convergence is quite good. For comparison, we present also a 
deviation of the theoretical ionization potential J of atoms He 4 from the experimental value, 
jexp ^ i983io.8cm _1 , by the data from [THl EEl EE] ■ It is seen that the convergence for the 
energy is faster, which is natural. Indeed, at the determination of namely the energy is 
minimized. However, the quantities (ff) and {ff) are determined also with a high accuracy. 
In work [26J, \l/o was calculated by means of the expansion in Laguerre polynomials. Up to 
1078 terms were taken into account, and the determined value of (ff) (see Table) coincides 
with accuracy to 4 decimal points with the value obtained by us for the 80-parameter ^/ - 
This means that the above-determined (ff), (ff), and the upper and lower limits of D-j for 
this \l/o have error less than 0.1%. The modern more exact methods of calculations of are 
available (see, e.g., (27j El]), but the too high accuracy is not required for our task, because 
it will not increase the accuracy of D 7 , since the large error of D-j remains unremovable due 
to the "smearing" of aE. 

Since the main contribution to the DM is given by the states with large n, for which 
aE w AE ioni D 7 should be closer to the lower limit. Therefore, we define D 7 in the table as 

D 7 w D r 7 nin + (D 7 nax - £>™ n ) /4. (36) 

In view of the result for the 80-parameter we finally have 

D 7 ^D 7 ^ 25.2 ± 2. (37) 

Solution ( |26l ). ( 1371 ) is in agreement with (jTJ), though we have obtained a somewhat larger 
value for D 7 . We note that the formula different from (1271) for D 7 was derived in [HI HE] 
with the use of dispersion sums. The most exact value, D 7 = 18.4, was obtained there on 
the basis of the 20-parameter \l/ taken from [23]. In this case, certain approximations were 
used in calculations, but the error arising due to these approximations was not evaluated. 
For this and more exact functions \l/o, we have calculated a more larger value of D 7 (l37l) . The 
difference of values of D 7 in ( 1371 ) and (pQ) is probably related to the approximations made in 
both calculations. We note that the above-proposed calculation is much simpler than that in 

pl HE]. 

3 Conclusions 

The above-considered (and earlier in [HI [IE]) tidal mechanism of the polarization of atoms 
is, in our opinion, a key to the comprehension of electric properties of superfluid He 4 and 
other fluids of nonpolar atoms, though we cannot yet exclude the contribution of other 
mechanisms (for example, the gradient thermoemf arising in an electrode due to the gradient 
of temperature). At the present time, the modeling of the near-surface polarization of helium 



near the electrode seems to be the most promising for the explanation of observed electric 
signals [HE]- And we do not exclude the possibility that the researchers of this field will meet 
a number of surprises. 
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